QPOs: Einstein's gravity non-linear resonances 
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Summary. There is strong evidence that the observed kHz Quasi Periodic Oscil- 
lations (QPOs) in the X-ray flux of neutron star and black hole sources in LMXRBs 
are linked to Einstein's General Relativity. Abramowicz&Kluzniak (2001) suggested 
a non-linear resonance model to explain the QPOs origin: here we summarize their 
idea and the development of a mathematical toy-model which begins to throw light 
on the nature of Einstein's gravity non-linear oscillations. 



1 Introduction 

QPOs are highly coherent Lorentzian peaks observed in the X-ray power 
spectra of compact objects in LMXRBs. They cover a wide range of fre- 
quencies (from mHz to kHz) and they show up alone, in pairs or in triplets. 
In the present review we focus on kHz QPOs which occur in pairs (the so 
called twin peak QPOs). Many models have been proposed in order to un- 
derstand kHz QPOs: some of them involve orbital motions (e.g. [14], [8], but 
there are also models that are based on accretion disk oscillations (e.g. [17], 
[7], [13], [9]). Here we present semi-analytical results, connected in particular 
with the Kluzniak&Abramowicz model. 



2 QPOs and General Relativity 

High frequency (kHz) QPOs lie in the range of orbital frequencies of circular 
geodesies just few Schwarzschild radii outside the central source. Moreover 
the frequencies scale with 1/M, where M is the mass of the central object 
[10]. These two facts support the strong gravity orbital oscillations models. 
Consider a test particle rotating around a compact source: the radial epicyclic 
frequency of planar motion and the vertical epyciclic frequency of nearly off- 
plane motion are respectively defined as 
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where U eff = g u + Ig^ + Pg^ is the effective potential and {r,6,4>) 
are spherical coordinates. The eigenfrcqucncies depend only on the metric 
of the system, hence on strong gravity itself. In Newtonian gravity there is 
degeneracy between these eigenfrequencies and the Keplerian frequency (all 
three frequencies are equal) while in General Relativity this degeneracy is 
broken and as a consequence two or three different characteristic frequencies 
are present, opening the possibility of internal resonances . As a consequence, 
while Newtonian orbits are all close, in GR they do not close after one loop 
(this is for the same reason as the well-known advance of the perihelion of 
Mercury). 

For a spherically symmetric gravitating fluid body (a better model of the 
accretion disk) these are the frequencies at which the center of mass (initially 
on a circular geodesic) oscillates. 

3 Kluzniak- Abramowicz resonance model 

QPOs often occur in pairs, and the centroid frequencies of these pairs are in 
rational ratio (e.g. [15]): these features suggested that high frequency QPOs 
are a phenomenon due to non-linear resonance, and that there may be an 
analogy between the radial and vertical oscillations in a Shakura-Sunyaev disk 
and the motion of a pendulum with oscillating point of suspension ([4]). Since 
in GR ui r < uj z the first allowed resonance would appear for ui z : uj r = 3 : 2 
and it would be the strongest. 

In all four microquasars which exhibit double peaks, the ratio of the two 
frequencies is 3 : 2, as well as in many neutron star sources. Moreover com- 
binations of frequencies and subharmonics have been detected: these are all 
signatures of non-linear resonance and they confirm the validity of the model. 

3.1 Toy model 

A mathematical toy-model for the Kluzniak- Abramowicz resonance idea for 
QPOs was recently developed [1]. It describes the QPOs phenomenon in terms 
of two coupled non-linear forced oscillators, 

Sr + uj 2 r Sr = F(Sr, Sz, Sr, Sz) + C cos (u t) + N r (t), (2) 

Sz + lo 2 z Sz = G(5r, Sz, Sr, Sz) + Dcos (u t) + N z (t), (3) 

where F and G are polynomes of second or higher degree (obtained in 
terms of expansion of the deviations from a Keplerian flow). The cos (uiot) 
terms represent an external forcing: they are mostly important in the case of 
NS, where wo can be the NS spin frequency. N r and N z describe the stochastic 
noise due to the Magneto-Rotational Instability . 

In absence of turbulence, the first finding, which can be derived by using 
any perturbative method, is that in the approximate solution there are terms 
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with the denominators in the form ruv r — mu> z (with n and m being integers). 
These n and m cannot take any value, but they depend on the symmetry 
of the metric and of the perturbation: for a plane symmetric configuration 
one can demonstrate (e.g. [12], [6]) that m = 2p (p integer). Due to this, the 
regions where nuj r = 2pu z are candidate regions of internal resonance , and 
the strongest one would be for lo z : lo t — 3 : 2 (a similar reasoning can be 
done for the forcing term). 



3.2 Correlation and anticorrelation 



A well known property of weakly non-linear oscillators is the fact that the 
frequencies of oscillation (u>* and w*) depend on the amplitudes, even if their 
value remains close to that of the eigenfrequencies ( u) r and uj z ). Perturbative 
methods help us again (e.g. [11]): indeed one can write 



W ; - uj z + e 2 ^ + e 3 ^ + 0(e 4 ) (4) 

and find the frequencies corrections {u>l ) by constraining the solution to 
be stationary at a given timescale (of the order of e~ J ). For small perturba- 
tions one would get that to* z — Au>* + B. In this way one can qualitatively 
explain the observed linear correlation between the twin frequencies in NS 
sources: the example of the NS source Sco X-l was studied by [3] and[12] (see 
Fig.l). 
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Fig. 1. The dotted line is the least- »W 

squares best-fit to the data points (the Fig _ 2 The anticorrelation between 

observed kHz QPOs frequencies in Sco shifts and glopes {A) The pointg 

X-l); the thin solid line corresponds to , , . , n , 

" correspond to the mdividuaf Z and 

a slope of 3 : 2 (for reference) . The atoU sourceg ^ [2] The begt 

thick solid segment is the analytic ap- fit ^ goeg through the pQint ^ L5) 
proximation, in which the frequencies (Courtegy of Qabriel T6r6k) 
are scaled for comparison with obser- 
vations. 
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A direct consequence of the internal resonance is the fact that for different 
sources the coefficients A and B should be nearly linearly anticorrelated ([1]): 
A = u) z /u) r — B/f (M), where fo{M) is a function which depends on the mass 
of the central object. Note that at higher orders the relation is expected to 
slightly deviate from linearity: anyway up to now this feature fits very well 
the available data (see Fig. 2). 

3.3 The effect of turbulence 

Accretion disks are characterized by a huge number of degrees of freedom. 
The turbulent processes can be assumed to have a stochastic nature. In par- 
ticular, we have investigated a simplified model for the Kluzniak- Abramowicz 
nonlinear theory and showed that a small noise in the vertical direction can 
trigger coupled epicyclic oscillations (see Fig. 3 and 4). On the other hand too 
much noise would disrupt the quasi-periodic motion [16]. This is similar to 
the stochastically excited p-modes in the Sun, and it may help in estimating 
the strength and the nature of the turbulence itself. 
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Fig. 3. Low turbulence: power spec- 
tra (upper part) and phase diagrams 
(lower part) for r and z. The displace- 
ments are in units of r g , the frequen- 
cies are scaled to kHz (e.g. assuming a 
central mass M of 2M ). It does not 
differ much from the behavior in ab- 
sence of turbulence. 




Fig. 4. The same as in the previous 
figure. In this case however the tur- 
bulence is strong enough to feed the 
resonance: as a consequence the ampli- 
tudes of oscillation are much greater. 



QPOs: Einstein's gravity non-linear resonances 5 

4 Conclusions 



Non-linear parametric resonances occur everywhere in Nature: together with 
GR they could explain the mechanism at the basis of kHz QPOs. In this 
way the mass and the angular momentum (e.g. [5]) of the central compact 
object could be precisely measured, but most of all Einstein's strong gravity 
could be proved: a good motivation to keep on investigating on this puzzling 
phenomenon. 

References 

1. M.A Abramowicz (editor), "Nordita Workdays on QPOs" :Astron.Nachr., 326, 
9 (2005) 

2. M.A Abramowicz et al. in: S. Hledfk, Z. Stuchlfk (eds.), Proceedings of RAG- 
tvme 6/7 , Silesian University Opava (2005) 

3. M.A. Abramowicz,V. Karas, W. Kluzniak, W.H. Lee, P. Rebusco: PASJ 55, 
467 (2003) 

4. M.A.Abramowicz, W.Kluzniak: A& A, 374, L19 (2001) 

5. B. Aschenbach: A & A, 425, 1075 (2004) 

6. J. Horak in: S. Hledfk, Z. Stuchlfk (eds.), Processes in the Vicinity of Black 
Holes and Neutron Stars, Silesian University Opava (2004) 

7. S. Kato:PASJ, 53,L37 (2001) 

8. F.K.Lamb, M.C. Miller: astro-ph/0308179 (2003) 

9. L.-X Li,R. Narayan: ApJ, 601, L414 (2004) 

10. J.E. Mc Clintock, R.A. Remillard in Compact Stellar X-ray Sources, ed. 
W.H.G. Lewin& M. van der Klis (Cambridge: Cam. Univ. Press), astro- 
ph/0306213 (2004) 

11. DA. Mook, A.H. Nayfeh: Nonlinear oscillations, John Wiley& Sons, New York 
(1976) 

12. P.Rebusco: PASJ, 56, 553 (2004) 

13. L.Rezzolla, S'i Yoshida,T.J. Maccarone, O. Zanotti: MNRAS, 344, L37 (2003) 

14. L.Stella, M. Vietri: ApJL, 492, L59 (1998) 

15. T.E.StrohmayenApJL, 552, L49 (2001) 

16. R.Vio, P. Rebusco, P.Andrcani,R.V.Overgaard:A&A (submitted 2005) 

17. R.V. Wagoner,A.S. Silbergleit, M. Ortega-Rodnguez:ApJ, 559, L25 (2001) 



